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CROSS REFERENCES 

[0001] The present application is co-pending with and claims priority benefit of U.S. 
provisional application entitled "Geolocation of Mobile Appliances" S.N. 60/418342 and filed 
on October 16, 2002, the entirety of which is incorporated herein by reference; 

[0002] The present application is related to and concurrently filed with applications titled 
"A NETWORK OVERLAY GEO-LOCATION SYSTEM WITH SMART ANTENNAS AND 
METHOD OF OPERATION" SN #, "WIRELESS COMMUNICATION NETWORK 
MEASUREMENT DATA COLLECTION USING INFRASTRUCTURE OVERLAY-BASED 
HANDSET LOCATION SYSTEMS" SN #, "NETWORK OVERLAY LOCATION SYSTEM 
AND METHOD FOR AIR INTERFACE WITH FREQUENCY HOPPING" SN #, "A SYSTEM 
AND METHOD FOR ENHANCING THE ACCURACY OF A LOCATION ESTIMATE", 
SN#, and "SYSTEM AND METHOD FOR OPERATING A NETWORK OVERLAY GEO- 
LOCATION SYSTEM WITH REPEATERS" SN #, filed October 16, 2003, the entirety of each 
of these applications is incorporated herein by reference. 

BACKGROUND 

[0003] Digital signals are often filtered using a pulse shaping filter prior to transmission. 
This is typically done to contain the signal bandwidth and minimize intersymbol interference 
between signal components corresponding to different digital symbols. This is shown in 
Figure 1, where symbols from a digital constellation corresponding to the information being 
modulated are passed though a pulse shaping filter. In Figure 1, the digital data to be transmitted 
is mapped into a complex signal constellation in block 101. For example, the complex signal 
constellation used may be an M-ary QAM constellation; however other constellations are also 
used. The mapped constellation undergoes pulse shaping in a filter as shown in block 102. 
Several methods known in the art can be employed for pulse shaping. The filtered constellation 
signal is converted to a radio frequency, represented as block 103, for transmission over the 
ether. 
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f 0004J An artifact of this signal generation provides an unintended feature that has been 
the focus of recent investigation. If the digital signal is passed through a multi-path channel, the 
channel output is an aggregate of delayed, possibly faded and phase shifted replicas of the 
original digital signal. In practice, this occurs if a multiplicity of reflections of the transmitted 
signal are contained in the received signal. The delays can either be absolute, if the time of 
arrival of the direct path signal is known, or could be relative delays between the multi-path 
components. If these multi-path signals are received at the antenna array, the received signal can 
be mathematically formulated as a space-time signal. 

[0005] When the characteristics of the pulse shaping filter and the antenna array are 
known, a theory of signal processing can be applied to estimate the multi-path delays of the 
signal components and their particular directions of arrival. This signal processing analysis is 
referred to as space-time processing. Space-time processing is a group of techniques that may be 
applied to resolve the received space-time signals into a sum of faded space-time signals. Each 
of these space time signals corresponds to the particular angle of arrival and time delay of one of 
the original multi-path signal components. 

[0006] It is advantageous to develop a mathematical description of the prior art technique 
to convey the manner in which multi-path delays and angle of arrival ("AO A") are currently 
calculated to fiilly appreciate the distinctness of the to be disclosed subject matter. The prior art 
method is illustrated in Figure 2. 

[0007] A column vector r k denotes the received signal at antenna £ of an antenna array 
with m antennas, where k = l,2,...m. An impulse response h k of the multi-path channel is 
derived from r k , represented in block 201. The derivation of the column vector h k can be 
achieved by various methods and implemented with signal processors through software and/or 
hardware. 

[0008] If the source data associated with this received block is known, a simple means of 
extracting h k is via the delay matrix corresponding to this source data. The delay matrix Z is 
formed by stacking symbol shifted copies of the source data in rows to a depth that defines the 
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extent of the desired impulse response and truncating its longer dimension to match the length of 
r k . An estimate of h k is given by: h k = (ZZ**)' 1 Zr k . 

[0009] Alternate means for estimating the impulse response may provide better or worse 
estimates, depending on the particular modulation format of the data, the block length, the fading 
characteristics of the multi-path channel and possibly other parameters. Some of these other 
methods are blind to the actual data transmitted, using properties of either the signal modulation 
and/or of the channel instead. 

[0010] Having estimated the impulse response of the multi-path channel from the source 
to antenna k of the array, a vectorized space-time impulse response over the entire array is 
formed, in block 202, by stacking the individual impulse response estimate hk into a long column 

vector /, given by: I = 



[0011] Theoretically, / can be expressed as / = ^J, + N 9 where 7 indexes the 

individual space-time impulse responses, i = 1,2,..., n, of the distinct multi-path components and 
N is a noise vector. 

[0012] Any particular /, is of the form /, = ftrjaiOJ ® s(r,) in which ft denotes the 
fade multiplier for the signal block, and 77 denotes the signal amplitude at the transmitter. a(0 t ) 
denotes the antenna response corresponding to a signal arriving from angle 0. , ® denotes the 
Khatri-Rao product, and g{r t ) denotes the pulse shaping waveform delayed by r ( and sampled. 
This formation of an outer product and aggregate in the covariance matrix, is represented in 
block 203. The formation of the covariance matrix can be implemented with signal processors or 
other computer processors through software and/or hardware devices. 

[0013] The vectors a and g can be expressed as: 



K • 
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, and 



g(-(/"l)7;-r t .) 



[0014] In the equation for g(r t ) , / denotes the sampling depth of the pulse shaping 
function and T s is the sampling time. 

[0015] Given this formation of the space-time impulse response, when the number of 
multi-path components is smaller than the dimension of the symbol sampled impulse response 
vector/, it is possible to estimate the multi-path delays r. and the multi-path arrival angles 0. . 

[0016] The prior art approach to estimating the delays and arrival angles relies on an 
explicit knowledge of the aggregate of all vectors: 



for all angles 6 X . This aggregate is termed the array manifold, A. It is 

K,(*.)J 

assumed that the pulse shaping function at the transmitter is known at the receiver. Denoting the 
aggregate of all vectors as 



g(-lT,-T t ) 

gt-V-iyr.-Tj 



for all values of x i as the delay manifold, g, then the quantity K = A® g represents the space- 
time manifold. 
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[0017] The observation that 7, is contained in K leads to a primary objective of space- 
time processing: searching the manifold AT for weighted linear combinations of vectors // such 
that a best fit to the observed space-time impulse response I is generated as shown in block 204. 
A variety of techniques may be applied for this purpose, such as Multiple Signal Classification 
(MUSIC), The Method of Alternating Projections (APM), etc, which can be implemented 
through software and/or hardware. Other mathematical descriptions for jointly estimating the 
angle of arrival ("AOA") and time delays can be found in Ziskind, I., Wax, M., "Maximum 
likelihood localization of multiple sources by alternating projection", IEEE Trans. Acoust, 
Speech, Signal Process, vol. 36, no. 2 (Oct. 1988), 1553-1560; Van Der Veen, M, Papadias, 
C.B., Paulraj, A.J., "Joint angle and delay estimation" IEEE Communications Letters vol. 1-1 
(Jan. 1997), 12-14; Schmidt, R.O. "Multiple emitter location and signal parameter estimation" 
Proc. RADC Spectrum Estimation Workshop, (Mar. 1999), 243-258; Young-Fang Chen, 
Michael D. Zoltowski "Joint Angle and Delay estimation ofDS-CDMA communication systems 
with Application to Reduced Dimension Space-time 2D Rake Receivers", IEEE Transactions on 
Signal Processing; Paulraj, A.J., Papadias, C.B., "Space-Time Signal Processing for Wireless 
Communications", IEEE Signal Processing Magazine, vol. 1 1 (Nov. 1997), 49-83; Paulraj, A. J., 
Papadias, C.B., "Space-Time Signal Processing for Wireless Communications: A Survey" 
Information System Laboratory, Stanford University; and Haardt, Brunner and Nossek "Joint 
Estimation of 2-D Arrival Angles, Propagation Delays, and Doppler Frequencies in Wireless 
Communications"; all of which are incorporated herein by reference. 

[0018] An object of the disclosed subject matter is to obviate the deficiencies of the prior 
art by removing the dependency of the time delay estimates from the spatial and gain 
characteristics, of an antenna array thus allowing multi-path delay estimates to be obtained for 
any generic antenna array. This object is achieved by recasting the array manifold in a spatially 
blind manner so as to be independent of the array characteristics. 

[0019] It is another object of the disclosed subject matter to present an improved method 
for estimating the multi-path delays in a signal received at any k array element. The method 
includes estimating an impulse response at each k antenna, generating a space-time impulse 
response, and forming a covariance matrix and resolving the covariance matrix with a known 
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antenna array manifold. Additionally, a novel improvement to known methods includes the step 
of resolving the covariance matrix with a fictitious antenna array manifold. 

[0020] It is still another object of the disclosed subject matter to present a novel method 
for estimating the multi-path delays in a signal using a spatially blind antenna array. The method 
includes generating an impulse response A* for each antenna k in the antenna array and 
determining a vectorized space-time impulse response / over the antenna array. The method 
further includes creating a covariance matrix C, a fictitious manifold^/, where A f is spatially 
blind and independent of the array characteristics, and then resolving the covariance matrix C 
with the fictitious manifold A f \o estimate the multi-path delays r, in a manner independent of 
the array characteristics. 

[0021] It is yet another object of the disclosed subject matter to present a method of 
estimating the multi-path delays of a sequence of / blocks of a signal received at an antenna array 
comprising k antenna elements independently of the spatial array characteristics of the antenna 
array. The method includes deriving a channel impulse response estimates Ay,* for each block j at 
each antenna k and determining a vectorized aggregate space-time impulse response / for each 
blocky. The method includes the steps of forming an estimated covariance matrix for the 
sequence of / blocks, forming an array manifold Af void of spatial information; and then 
resolving the covariance matrix with the fictitious array manifold Af to determine the multi-path 
delays . 

[0022] It is also an object of the disclosed subject matter to present a novel system for 
estimating the multi-path delays in a signal using a spatially blind antenna array. The system 
includes an antenna array, a means for generating an impulse response ht t a means determining a 
vectorized space-time impulse response / and a means for creating a covariance matrix C. The 
system also includes a means for creating a fictitious manifold Af, wherein Af is spatially blind 
and independent of the array characteristics; and a means for resolving the covariance matrix C 
with the fictitious manifold ,4/ to estimate the multi-path delays r i independent of the array 
characteristics. 
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[0023] These objects and other advantages of the disclosed subject matter will be readily 
apparent to one skilled in the art to which the disclosure pertains from a perusal or the claims, the 
appended drawings, and the following detailed description of the preferred embodiments. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0024] Figure 1 is a representative flow chart of a portion of a prior art system and/or 
method for transmitting a digital signal. 

[0025] Figure 2 is a representative flow chart of a prior art receiving system and/or 
method for space-time processing a received signal. 

[0026] Figure 3 is a representative flow chart for spatially blind space-time processing a 
received signal according to a system/method of the present disclosure. 

[0027] Figure 4 is a representative flow chart for estimating delays of multipath 
components of a received signal according to an embodiment of a system/method of the present 
disclosure 

DETAILED DESCRIPTION 

[0028] The present disclosure presents a system and method for a technique for 
determining the individual times of arrival of multipath signal components in a received signal. 
The delays in the multipath components are determined by creating a fictitious array manifold 
A f in which the directional knowledge of the array is absent. As discussed earlier, the vector a 

is produced with the knowledge of all the vectors for all angles 0 i , is represented as: 
«a(*i) 

[0029] To remove the dependence on $, , each dimension of the array takes all possible 
complex numbered values. The independent aggregate vector is then given by: 



a(0,) = 
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, where each a k is contained in the complex number space C, so that the 



fictitious array manifold A f is the space C m . This formulation allows space-time processing to 

proceed blind to the specific characteristics of the antenna array such as element spacing, 
element gain, etc. 

[0030] As discussed above, the signal processing flow for the established prior art 
technique is shown in Figure 2. This technique requires knowledge of the antenna array 
manifold A, i.e., one must have explicit knowledge of the antenna response corresponding to a 
signal arriving at each possible angle for each element in the antenna array. As discussed earlier, 
it is not always possible to easily obtain the antenna array manifold A as calibration of the 
antenna array is a tedious, time-consuming, and expensive task. The prior art technique allows 
for determination of an estimation of multipath delays and the particular angles of arrival of the 
multipath signals. 

[0031] The signal processing flow for the above-discussed subject matter, on the other 
hand, employs a fictitious array manifold as shown in Figure 3. An estimate of the impulse 
response is determined at each antenna element as shown in block 301. From the impulse 
response, a space-time impulse response is generated in block 302 and an outer product and 
aggregate in the covariance matrix is formed in block 303 in a manner similar to that as 
discussed for the prior art. In block 304 the covariance matrix is resolved with the generalized 
fictitious antenna array manifold A/ defined above rather than the known antenna array manifold 
A as shown in Figure 2. While the novel techniques described herein does not allow for 
resolution of the particular angles of arrival of the multipath signals, it does allow for 
determining an estimation of the time delays of the multipath components given an antenna array 
with arbitrary and unknown properties. Therefore, the novel techniques herein described allow 
for the use of antennas separated by arbitrary distances and with arbitrary gain characteristics 
while still obtaining important information: the time delays of multipath signal components. 



9 



WO 2004/036525 



PCT/US2003/032585 



[0032] To facilitate understanding of the novel techniques described herein, an example 
follows for the case of a signal received at a two antenna array in which the established prior art 
space-time techniques derive delay estimates on a multi-path signal, followed by estimates 
developed by use of a spatially blind antenna array according to an embodiment of the present 
disclosure. 

10033] A sequence of blocks of signal received concurrent in time at the two antennas are 
denoted by rjk where k = 1, 2 and /= 1 ,2, ... J, the index j counts the blocks and of course the index 
k references the antennas, two in this case. It is, of course, to be understood by those of skill in 
the art that the present discussion is exemplary only and that the techniques described herein are 
in no way limited to antenna arrays composed of only two antenna elements. The disclosed 
techniques are applicable to antenna arrays of an arbitrary number of antenna elements. If the 
source data associated with each block is known, possibly from demodulating the received signal 
at either or both of the antennas, or by demodulation at some other receiver, the delay matrix can 
be developed as illustrated in the previous section and this can be used to derive channel impulse 

response estimates for each block at each antenna. Thus h JJt = (ZjZj** )" ! Zj r j.k » w here h jk is the 
impulse response estimate for block j and antenna k, and Zj is the delay matrix for block j. 

[0034] The vectorized aggregate space-time impulse response for block j is given by 

/ , = ^ JtX , from which an estimated covariance matrix for the sequence of blocks can 
j 

then be formed as C=^IjI H . 

M 

[0035] Most known methods for extracting an arrival angle do so by operations on C. 
Other methods that operate on the sequence {Ij) differently can also be used. 
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[0036] All established techniques for estimating the path delays and arrival angles in the 
space-time context require knowledge of the array manifold A which is the aggregate of all 
possible vectors: 

in which a k {$) is the complex gain of antenna k for a signal arriving from angle 0 . Explicit 
knowledge of antenna gain for each antenna is required. In addition, the pulse shaping function 
and thus the delay manifold gas detailed in the previous section are also needed. 

[0037J Most established prior art techniques estimate the delays and arrival angles of the 
multi-path signals by decomposition of the covariance matrix into a sum of estimated impulse 
responses. This is done using the space-time manifold: 

K=A<2> g. 

[0038] The estimated covariance matrix is resolved into estimates of the impulse 
responses associated with each multi-path component. Each estimated impulse response is 
drawn from K and therefore points directly to a particular delay and arrival angle. 

[0039} In accordance with an embodiment of the disclosed subject matter, an example 
using a spatially blind estimation techniques for the time delay estimates is illustrated below for 
the same scenarios as described above. 

[0040] To obtain the delay estimate without knowledge of the antenna array or an 
explicit knowledge of the antenna gains , a fictitious array manifold is used, A/, given by the 

aggregate of all vectors, a = | 1 where a, and a 2 range over the set of complex numbers. It is 

worthy to note that while the delay estimates can be determined for multi-path signals 
independently of the antenna characteristics, angle of arrival cannot be determined in the same 
manner. This is not a problem if all that is of concern is the relative time delays of the multipath 
signals and the angles of arrival of those multipath signals is not of consequence. 
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[0041] The process estimates the delays without knowledge of the actual array manifold 
A by using instead the fictitious manifold A/, so that A/ replaces A in the space-time manifold 



space of complex numbers and r covers the expected range of multi-path delays. 

[0042] The estimated covariance matrix is resolved in order to estimate the impulse 
responses associated with each multi-path component. Each estimated impulse response is 

drawn from K f and maps directly to a particular r and complex vector \a x , a 2 ] . The signal 

processing flow for this embodiment is shown in Figure 4. An estimate of the impulse response 
is determined at each antenna element as shown in block 401 . From the impulse response, a 
space-time impulse response in generated in block 402 and an outer product and aggregate in the 
covariance matrix is formed in block 403. In block 404, the covariance matrix is resolved with 
the specific fictitious antenna array manifold A/. 

[0043] The earliest arrival time is of particular interest in many applications such as 
cellular geo-location. This technique, termed BLind Space Time ("BLST") permits the 
estimation of relative delays in multi-path signal components using arbitrary configurations of 
antennas in an array. This, for example, allows one to apply the method to antennas that are 
separated by arbitrary distances and with arbitrary gain characteristics. Existing antenna 
configurations in the cellular wireless network can be used without the need for calibration of the 
array. Calibration is a tedious task that needs to be updated often as a known array is required 
for signal processing according to the prior art. Thus the use of Blind-Space Time for estimating 
multi-path delays allows for the application to any arbitrary antenna array without requiring 
calibration. 

[0044] While preferred embodiments of the present inventive system and method have 
been described, it is to be understood that the embodiments described are illustrative only and 
that the scope of the embodiments of the present inventive system and method is to be defined 
solely by the appended claims when accorded a fUll range of equivalence, many variations and 
modifications naturally occurring to those of skill in the art from a perusal hereof. 



K = A® g . This then becomes the aggregate of vectors, K = 




as a, and a 2 cover the 
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